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Abstract 

Cloud base height profiles measured with laser ceilometer are studied using multi- 
fractal approach. The irregular structure of the signals is a benchmark for nonlinear 
dynamical processes. A hierarchy of generalized dimensions determines the intermit- 
tency of the signal. The multi-affine properties are described by the h{-y) function. 
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1 Introduction 



Cloud base height profiles (CBHP) are known to have highly fluctuating and 
irregular structure. The dynamics of CBHP evolution is determined by a va- 
riety of processes in the atmosphere [1]. This irregular structure is a bench- 
mark for nonlinear dynamical processes. We apply the multifractal approach 
to study the scaling properties of the CBHP. In this report we present an 
analysis of the multi-affine properties of a CBHP signal y{t) (Fig.l) measured 
with a ground-based laser ceilometerQ having a temporal resolution of 30 sec- 
onds, taken on September 23-25, 1997 at the Southern Great Plains (SGP) 
(Oklahoma, USA) site of the Atmospheric Radiation Mesurement Program of 
the Department of Energy. 

^ http://www.arm.gov 
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Fig. 1. Cloud base height profile data measured at SGP site on Sept. 23-25, 1997. 
Time resolution : 30 sec; 7251 data points. 

2 Multi-afRnity and Intermittency 

First, we tested the scaling properties of the power spectral density S{f) 
of the CBHP signal and obtained (figure not shown) that S{f) with 
(5 — 1.46 ± 0.06 for frequencies lower than 1/15 min~^ and /3 for higher 
frequencies. 

The multi-affine properties of y{t) can be described by the so-called "q-th" 
order structure functions [2-4] Cq = {\y{ti+r) — y{ti)\'^) {i = 1) 2, . . . , — r), 
where the averages are taken over all possible pairs of points that are r — 
ti+r—ti apart from each other with r > 0. Assuming a power law dependence of 
the structure function, the H{q) spectrum is defined through [2] Cq{r) ~ r^^^^^ 
with g > 0. 

The intermittency of the signal can be studied through the so-called singular 
measure analysis. The first step that this technique requires is defining a basic 
measure e{l; I) as e{l; I) = \ Ay{l; l)\/ < Ay{l; /)>(/ = 0, 1, ... , A^-1), where 
Ay[l;l) = — y{ti) is the small-scale gradient field and < Ay{l;l) >= 

jfYliLo^ |A|/(1;/)|. It should be noted that we use spatial/temporal averages 
rather than ensemble averages, thus making an ergodicity assumption [5] as 
our only recourse in such an empirical data analysis. Next we define a series 
of ever more coarse-grained and ever shorter fields e{r;l) where < / < 
N — r and r = 1,2,4, ...,iV = 2"^. The average measure in the interval 
[/; l + r]is e{r; I) = ^ EjtV^ e{l] V) {I = 0, . . . , N - r) The scaling properties 
of the generating function are then searched for through ^^(t) =< e{r;iy >~ 
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q--K{q) ^ for g > 0. Thus the multi-fractal properties of the CBHP signal are 
expressed by two sets of scaling functions, H{q) describing the roughness of the 
signal and K[q) describing its intermittency. The intermittency of the signal 
can be also expressed through the generalized dimensions D{q) as introduced 
by Grassberger [6] and Hentschel and Procaccia, [7] D{q) — 1 — K{q)/{q — 1). 




Fig. 2. Hierarchy of generalized dimensions D{q) for the CBHP data in Fig.l. The 
straight line is drawn to enhance the value Qs at which D{q) function starts to 
deviate from a linear dependence. This deviation is related to sampling problems 
due to a single realization. The dashed line defines the monofractal case Di = 0.905. 



The multi-affinity of y{t) means that one should use different scaling factors 
H{q) = Hq in order to rescale such a signal. This also implies that local 
roughness exponents 7 exist [8] at different scales. The density of the points 
N^{t) that have the same roughness exponent is assumed [9] to scale over the 
time span r as N^{t) ~ T~h{'j). From Ref. [10] the following relations are 
found : 7(g) = d{qH{q))/dq and hip/^) = 1 + g7(g) — qH{q). In Fig. 3 the 
CBHP multi-affine properties are presented via the /i(7)-curve. 



3 Conclusions 



We have demonstrated the multi-affine structure of cloud base height profiles. 
Further work will be directed toward relating these statistical parameters to 
the dynamical properties of the clouds, an important step toward understand- 
ing, modeling and predicting their dynamical behaviour. 
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Fig. 3. The /i(7)-curve for the CBHP data in Fig.l 
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